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Abstract 

In this paper we compute the decay constant of the pseudo-scalar heavy-light mesons 
in the heavy quark effective theory framework of QCD sum rules. In our analysis 
we include the recently evaluated three-loop result of order for the heavy-light 
current correlator. The value of the bottom quark mass, which essentially limits 
the accuracy of the sum rules for B meson, is extracted from the nonrelativistic 
sum rules for T resonances in the next-to-next-to-leading approximation. We find 
stability of our result with respect to all types of corrections and the specific form 
of the sum rule which reduces the uncertainty. Our results /b = 206 it 20 MeV 
and fo = 195 it 20 MeV for the B and D meson decay constants are in impressive 
agreement with recent lattice calculations. 

PACS numbers: 12.38.Lg, 14.40.Nd, 14.40.Lb 



1 Introduction 



The decay constant of a pseudo-scalar meson with one heavy and one light quark con- 
stitutes a hadronic quantity which is of primary phenomenological interest. It covers the 
strength of the leptonic weak decays of B and D mesons and enters as an input quantity 
into the analysis of the nonleptonic B and D meson decays and the B-B mixing pro- 
cess. The latter is of special interest since it provides a direct source of the information 
on the Cabibbo-Kobayashi-Maskawa matrix elements involving the top quark. Still no 
experimental information of sufficient accuracy is available and the theoretical study of 
the decay constant is mandatory. 

The first quantitative evaluation of the B and D meson decay constants Jb and fn 
using radiative corrections was performed in [|l| where the QCD sum rules proposed 
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in 0, 3 have been used together with the perturbative two-loop results of 0{as) 0. 
More refined evaluations followed ||, ^ ||, H where it was realized that the accuracy 
of the decay constants is significantly limited by the uncertainty of the bottom quark 
|ll], |12| the heavy quark effective theory (HQET) |jT3|, 



mass. In IT 



has been used 

to resum the leading and next-to-leading logarithms of the heavy quark mass. After the 
renormalization group improvement the two-loop corrections were found to be huge which 
might be a signal of the importance of the higher order contributions and questions the 
reliable determination of the decay constants. 

The only method to compute hadronic matrix elements which is entirely based on 
the first principles of QCD are probably the lattice gauge theory simulations. For the B 
meson decay constant the value 



lat 
B 



200 ± 30 MeV . 



is cited in the review article [0 as an average over different calculations performed on 
the lattice. It agrees well with a recent evaluation of fs [0 where the result fs = 
204(8) (29) (-1-44) MeV is given. Here, the first error is statistical, the second due to the 
discretization and the third one includes the uncertainty from the lattice scale. The 
numbers obtained for fs from lattice calculations are in reasonable agreement with the 



ones obtained from sum rules. An averaged value for the latter can be found in |17 



B 



178 ±42 MeV. 



For the D meson decay constant the situation is similar. The lattice result |T5| , \18 

f^' = 225 ±30 MeV, 
again overshoots the sum-rule estimate 



s.r. 
D 



188 ± 48 MeV 



(2) 



(3) 



(4) 



given in Ref. |T^. In spite of the essential progress in the lattice calculations their uncer- 
tainty is still rather large and it is too early to rely solely on the results of this approach. 
On the other hand, the sum-rule analysis can be essentially improved to reach an accuracy 
which is comparable or even better than the current accuracy of the lattice calculations 
and provides an independent cross check of the latter. In this paper we perform the anal- 
ysis of the sum rules for the pseudo-scalar heavy-light meson decay constant and improve 
the existing calculations with respect to several points which we summarize in the follow- 
ing: 

(i) We include the three-loop perturbative corrections of 0(«^) which recently became 
available |T|, |g. 

(ii) The bottom quark mass, which constitutes a crucial input for /b, is determined using 
the approach suggested in [^, in the context of semileptonic decays of the bottom 
quark. The latter is also quite sensitive to the pole mass rrib- The basic idea is to replace 
rrih by the ratio of experimental and theoretical moments of the nonrelativistic sum rules 
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for T resonances. It is computed order by order up to the next-to-next-to-leading accu- 
racy depending on the perturbative input used for the evaluation of the sum rules for fs- 
In this context we would like to refer to |^ where also the T-resonance sum rules have 
been used, however, only up to order as- 

The charm quark pole mass mc, necessary for the calculation of J'd, is then extracted 
from the HQET relation between mj, and rric- 

(iii) For the evaluation of the decay constant Jb we use both Laplace and Hilbert sum 
rules. We present, for the first time, the explicit formulae of the latter in the framework 
of HQET. The comparison of the results obtained with both approaches provides us with 
an estimate of the intrinsic uncertainty of the method. 

We will show that these new ingredients improve significantly the reliability of the 
sum rules and finally the prediction for the decay constant, in particular for fs- 

This paper is organized as follows. In Sections |] and |^ we introduce the basic features 
of the sum rules and the HQET formalism, respectively. In Section |^ the renormalization 
group improved sum rules for the decay constant are given in the three-loop approxima- 
tion. The problem of the heavy quark mass determination is discussed in Section ^. In 
Section ^, we present the numerical analysis and Section ^ contains our conclusions. 

2 Sum rules 

The decay constant of a pseudo-scalar meson P consisting of a heavy {Q) and a light 
quark (g) is defined through the matrix element 

mj;\pip)) = if pp., (5) 

where j'^ = q'y^'y^Q is the axial-vector current. In this paper we focus on the B meson 
with the bottom quark as the heavy constituent but present also the analysis of D mesons 
with charm as heavy flavour. We neglect the SU (3) violating effects of the strange quark 
mass. The ratio of the decay constants of strange and non-strange mesons can be reli- 
ably computed both on the lattice and using QCD sum rules as an essential part of the 
uncertainties cancel [|r^, [T^, [T^ . 

In order to derive the QCD sum rules for the pseudo-scalar heavy-light meson decay 
constant fp one considers the correlator 

^\q') = zJdxe''^^0\Td%{x)d^fJmO), (6) 

which is related to the correlator of the pseudo-scalar currents 

IP>{q^) = i Jdxe'''''{0\Tf{x)f\0)\0) (7) 

by the equation d^j'^ = rngj^, where = i{mQ{fi) /mQ)qj^Q and rfiQ^fi) and mg are the 
MS and pole mass of the heavy quark Q, respectively. 

Following [^] the correlator \E''^(g^) is evaluated in two ways. In the Euclidean region 
where q^ — mn <^ it can be reliably computed in QCD because of the asymptotic 
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freedom. The correlator gets a perturbative contribution corresponding to the leading 
operator in the operator product expansion (OPE) of the two currents in Eq. (|^). Fur- 
thermore there are power suppressed terms from the vacuum expectation values of the 
higher dimension operators (so-called vacuum condensates) parameterizing the nonper- 
turbative long-distance effects. On the other hand, the correlator can be obtained from 
the dispersion integral over the physical states. In the standard analysis only the ground 
state meson is supposed to give a delta-function contribution to the spectral function. 
Assuming local quark-hadron duality the contribution from the higher resonances is mod- 
eled by the perturbative continuum starting at some threshold Sc which brings an intrinsic 
uncertainty to the approach. Finally one arrives at the following equation 



nt{<l') + KM) = + d3 L -J + subtractions , (8) 

^ ^ Mp — TT s — q-^ 



where Mp is the meson mass. \l/pj and \l/^pt are the perturbative and nonperturbative 
QCD contributions, respectively. The subtractions needed for the dispersion integral are 
not specified explicitly as they will drop out in the following. To suppress the contri- 
bution from higher resonances and to reduce the uncertainty one can perform a Borel 
transformation of Eq. (Rl) 



1 / d 



(ra — 1)! \ dq 



— q^ ,11. — * oo TT 



1 foo p-s/M^ 



where M is the Borel parameter, and arrive at so-called Laplace sum rules. Alternatively, 
it is possible to consider moments of Eq. (||) 




^«(g2) 



(10) 



which leads to Hilbert sum rules. To estimate the intrinsic uncertainty of the method 
we will follow both options. In the above sum rules the weight functions cut off the 
dispersion integral at a typical hadronic scale much less than the heavy quark mass so 
that it is saturated by the near threshold region where the heavy quark is nonrelativistic. 
Laplace sum rules are particularly relevant for the nonrelativistic HQET analysis because 
of the heavy quark mass independent exponential suppression of the relativistic momen- 
tum region. The situation is more tricky in the case of Hilbert sum rules as is discussed 
below. 



3 Heavy quark effective theory 

Systematic description of the heavy quark nonrelativistic dynamics and consistent sep- 
aration of the relativistic effects can be done within HQET. Let us discuss this issue in 
more detail. The perturbative heavy-light quark system involves two dynamical scales: 
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the hard scale given by the heavy quark mass and the soft scale given by the off-threshold 
energy 



ijj 



1 1 



By integrating out the relativistic hard modes with the off-shell momentum of order mq 
one arrives at HQET which includes u as dynamical scale. The effect of the hard modes 
is accumulated in the Wilson (matching) coefficients leading to an expansion in along 
with relativistic corrections and contributions from higher dimensional operators leading 
to an expansion in I/itlq. In the hadronic matrix elements the latter is converted to 
an expansion in the dimensionless parameter A/mg where A Mp — uiq describes the 
nonperturbative long-distance effects and remains finite as itlq oo. In the process of 
scale separation spurious divergences appear at the intermediate steps which result in 
the anomalous dimensions of the effective theory operators and lead to the corrections 
involving the large logarithms of the form ln(mQ/a)). These logarithmic corrections can 
be resummed by solving the effective theory renormalization group equations. 

For this purpose let us consider the effective theory realization of the axial-vector 
current. The corresponding connection between the QCD operator and its HQET coun- 
terpart is given by 



j; = a(mQ)C(mQ)j;(mQ) + 0(l/mQ) 



(12) 



The matching coefficients have been computed in ^ up to order a^. In the MS 
subtraction scheme they read 

Ca{MQ) 



TT 



683 Uir^ 11 ^ ^ 

In 2 C 3 

576 72 18 36^^ ' 



n (— — 

^ \ 288 ^ 36 , 



1 + 



^9 /a("')(mQ)~ 



864 



TT 



(13) 



where ni is the number of light flavours and ({3) = 1.2020569 ... is the Riemann (- 
function. The renormalization group equation which governs the evolution of j^(yu) is of 
the form 



d/i2 

where the anomalous dimension is known up to two loops |^ 



(14) 



7 
7o 



7o- 

1 
2 



a 



(ni) 



TT 



+ 71 



a 



TT 



127 7tt^ 5 

7l = 1 Tl] . 

' 144 108 72 



(15) 
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The solution of Eq. (0) reads 



where 



270 //9o 



TT 



(16) 



(17) 



and the first two coefficients of the /^-function are given by 

11 1 ^ 51 19 

Thus in the nonrelativistic region u <^ mg we have the following representation of the 
perturbative part of the correlator (^ 



Im 



{Ca{mQ)C{mQ) 



2 Xjfi) 



mglm 



npt(cD) 



+ Oil/mQ) 



(19) 



where the universal HQET current correlator npt(ti;) does not depend on itlq and the spin 
and parity of the currents. Its imaginary part is known up to the three-loop approximation 
of Oial) 



3uj 



~,2 



1 



/1657 97772^ 



17 Air^ 



9 



L,. 



+ 



TT 



99(15) 



+ 



V 72 



54 



Lc, + —Ll 



nil 
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(20) 



where L^j = ln(/i2/a;2). The uncertainty in the non-logarithmic three-loop terms results 



from the semi-numerical method used in Ref. pO[. For jj ^ u the HQET correlator 
does not include large logarithms. They are all contained in the factor X{^)/X{mQ) 
which sums up the leading and next-to- leading logarithms of the form a"ln"(mQ/cj) and 
<+Mn"(mQ/c:;). 

On the phenomenological side the HQET decay constant is defined through the matrix 
element of the HQET current 



{o\]M\Hv)) 



/J. } 



(21) 



where \P{v)) is the nonrelativistic, i.e. quantum mechanical, meson state with velocity 
V. It is connected to fp via the relation 



fpJMp = Ca{mQ)CimQ)fp{mQ) + Oil/mQ) 



(22) 



By using the renormalization group property of the HQET current it is convenient to 
introduce the renormahzation group invariant quantity 

fp = V^/p(/i), (23) 
which is a universal low-energy parameter of strong interactions. 



4 The heavy quark effective theory sum rules 

Now we are in the position to write down the renormalization group improved sum rules 
for the heavy-light pseudo-scalar meson decay constant. 



4.1 Laplace sum rules 

Let us start with the sum rules in the infinite heavy quark mass limit mq —>■ oo. Trans- 
forming Eq. (^ to HQET and neglecting the mass suppressed terms we obtain 

ml{qq)\ 



TC JO 



dze-nm[U^,izT)]-{qq)ifi) 



1 + 2- 



TT 



+ X{T)- 



4T2 



(24) 



where A = (M|, - m'^)/mQ, Uc = (sc - rn?Q)/mQ, ml = {qgsG'"'cr^i,q) / {qq) and rngT = 
M^. We keep the operators up to dimension five in the OPE and neglect the running of 
the quark-gluon operator. 

The mass suppressed contribution to the correlator can be found by subtracting the 
asymptotic HQET result from the full theory expression. In this way the mass suppressed 
contribution to the physical decay constant from Laplace sum rules is obtained as 



+ {qq){mQ) 



2a, T 



mliqq) ^ {a^G^'^G^,) 



where 



Im 



Im 



TT rriQ Jo l + zT/mq] 2TmQ ' 127rmQ 

2 X(/.) 



Im 



npt(cu) 



The one-loop expression for this function is given by 

3 uj 



Im[5(^)H^,((:;)] 



^vrmg 1 + u/mQ 



The two-loop approximation for Hpt(g2) is known in full QCD in analytical form 



determines the 0{as) part of Eq. (^) which reads 



Im[5(2)H^p,(^)] 



2tt'^ 1 + 00 /mq 



UJ 



+ F 



UJ 



(25) 
(26) 

(27) 
|. It 

(28) 



where 



F{x) = 2Li2(-x) + ln(x) ln(l + x) - ln(x) + ^— ^ln(x + 1) - 1 

= Q + ^ln(x))x2 + 0(a:3), (29) 

with Li2{z) being the dilogarithmic function. The perturbative mass suppressed contri- 
bution of O(a^) can be obtained from the numerical three-loop result of |]19|, which 
is available under the URL http://www-ttp.physik.uni-karlsruhe.de/Progdata/ 
ttpOO-25. Thus we get the final expression for the decay constant 

fl . (^r ^£^t^^£^ijf , (30) 



MpJ Ximg) Mf 

The first term of this equation includes the leading HQET contribution up to C(«^) 
and the resummed leading and next-to-leading logarithms of the heavy quark mass. The 
second term represents all heavy quark mass suppressed terms up to order a^. 



4.2 Hilbert sum rules 

The use of HQET for Hilbert sum rules is rather subtle as they do not have a proper 
infinite heavy quark mass limit. Indeed, the naive limit rriQ oo leads to the decay 
constant 

(/;)' = ±^Col - (qq) + 0{a,) , (31) 

which, in contrast to Eq. (0), does not contain a dynamical constraint on the parameter 
Uc- Therefore one cannot use them to study the HQET decay constant. Nevertheless it 
is possible to apply Hilbert sum rules for the calculation of the physical decay constant 
and, furthermore, use HQET for the analysis. Indeed, if we keep the factor 1/s"^^ = 
l/(mQ""^^(l + cj/mg)""'"^) unexpanded in I/ttiq, it is straightforward to obtain the proper 
scaling 

= ^A + OiasA/mg) , (32) 

where A is defined after Eq. (|2^) , from the ratio of two arbitrary moments. Note that 
Eq. ( P^ is obtained from the purely perturbative correlator and the quark condensate 
contribution is neglected. This allows for the nonrelativistic treatment of the heavy quark 
in ^''^(s). Moreover, the dispersion integral is saturated by the region uj < mq/n and thus 
for n > niQ/ujc the result is not sensitive to Uc in contrast to Eq. (0). In this way we 
obtain the renormalization group improved Hilbert sum rules of the following form 

fp = -/ IV— ^1 — vm "^Q^ "^Q ) TT^^"" + 

[ TT Jo (1 + zujclmQ)'^+^ \ ^ ' X{mQ) 
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3 7r \^ niQ Jo (l + ztOc/mQy 



n{ n + 1) mliqq) ^ (a^G^^^) 
8 rriQ 127rmQ 



5 The heavy quark masses 

Before turning to the numerical analysis we want to discuss the determination of the 
heavy quark mass which is an input parameter of the sum rules given above. The result 
for the physical decay constant (^) is rather sensitive to the heavy quark mass value. 
Therefore mq should be determined with a great accuracy in order to obtain a reasonable 
precision for fp. The best accuracy of the bottom quark mass determination is achieved 
within the heavy quarkonium sum rules The corresponding expression for the pole 
mass is given by the ratio 



.Mr 



1 

2n 



= I • (34) 



Here the dimensionful experimental moments 



= 9/ d.:^^, (35) 





are generated by the normalized cross section of e~^e~ annihilation i?'^^P(s) = a{e^e^ — > 
hadronSb5)/(T(e"''e~ fi^fi^). The dimensionless theoretical moments are defined as 
follows 

A^:f = UniAmlrfds'^^^^, (36) 

where the vector current correlator is defined through 

(gV-gV'^)n''(g^) = z/da;e^^^(0|T/(a;)/(0)|0), (37) 

with Jn = Irj^h. If n is large enough the experimental moments are saturated by the T 
resonance contributions which is known with high precision. For large n the dispersion 
integral in Eq. (^) is saturated by the region near threshold where the nonrelativistic ex- 
pansion in the heavy quark velocity is applicable and the correlator can be systematically 
computed within the effective theory of nonrelativistic QCD (NRQCD) [^. The com- 
plete result for the moments including the second order corrections in the strong coupling 



constant and heavy quark velocity is now available ||2^, |2^, ^ 



It is widely believed that due to the renormalon contributions the absolute value of 



the heavy quark pole mass obtained through Eq. (|3^) is divergent [g^, ^ |3T[. As a 



consequence the absolute value of the pole mass is plagued with an intrinsic uncertainty 
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of order Aqcd- On the other hand, rrib is not an observable and has no immediate physical 
meaning. Therefore it can safely be removed from relations between physical observables. 
Using this philosophy we replace nib in the sum rules for Jb by the fixed order expression 
of the right-hand side of Eq. (p^ . Equivalently, we determine the value to the pole mass 
according to Eq. ( |34D only in a given order of the perturbative expansion correlated to 
the order of the approximation for fs- A detailed discussion of the sum rules and the 



corresponding numerical results can be found in |22[. In particular, we use the next-to- 
leading order (NLO) result, rrib = 4.68 GeV, for the calculation of J'b to order and the 
next-to-next-to-leading order (NNLO) result, nib = 4.79 GeV, for the C(a^) analysis of 
fs- The accuracy of the numerical value for the fixed order approximation for nib is no 
longer restricted by Aqcd since it is not related to the divergence of the series in Eq. (|3^ . 
It is mainly due to the dependence of the theoretical moments on the normalization scale 
of as and due to the dependence on n which altogether amounts to ±60 MeV for the 



NNLO result [El, 12^. Note that the direct order-by-order matching of the results for mb 



and /b is not obvious due to different kinds of resummation adopted in the sum rules: the 
study of the T-resonance sum rules requires the resummation of the singular Coulomb 
terms while the above analysis of the i?-meson sum rules involves the resummation of the 
heavy quark mass logarithms. Our matching of the perturbative series is based on the fact 
that in a given approximation both the T-resonance sum rules and the S-meson sum rules 
include the perturbative corrections of the same order in ag to the threshold behaviour 
of the heavy-heavy and heavy-light quark current correlators, respectively. Performing 
the analysis in the described way we expect that the large perturbative corrections to 
nib cancel in the complete expression so that the final series for Jb in terms of physical 
moments ( pS]) is convergent. This approach turned out to be very efficient for the analysis 
of the bottom quark semileptonic decay width We will show that the method 

works in the case of the sum rules for fs as well. 

For a given value of the bottom quark mass the charm quark mass can be obtained 
from the HQET constraint of the form 

mb-m, + 0{l/mb,c) = Mb - Md , (38) 

which results in mc = 1.37 GeV for the NNLO value of m^. Due to the cancellation 
between the different terms of order l/m^ c (see, for example, [^) this numerical value 



is valid with 0{l/mb^c) accuracy. The use of the relation (38) brings an additional uncer 



tainty to nic so that the total uncertainty can be roughly estimated as ±100 MeV. 

6 Numerical analysis 

In this section we present the numerical analysis of the sum rules. We adopt the same 
input values for the vacuum condensates as in ||T7| : 

{qq) (1 GeV) = - (225(25) MeV)^ , 
(asG^'^G^,) = 0.04(2) GeV\ 

mlilGeV) = 0.8(2) GeV^ (39) 
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The strong coupling constant is evaluated with four active flavours using two-loop ac- 
curacy and A*^'^'' = 296 MeV. This value corresponds to af'\mh) = 0.210 obtained from 
af'\Mz) = 0.1185 using the four-loop renormalization group evolution [Q. For the 
meson masses we use Mb = 5.2793(7) GeV and Md = 1.8641(10) GeV, respectively Q. 



6.1 The decay constant within heavy quark effective theory 

Let us start with the analysis of the limit tuq oo. The general philosophy for the 
determination of the decay constant from Laplace sum rules is as follows one has to 
optimize the upper bound of the duahty interval, Uc, in such a way that the value of fp 
as computed from Eq. (24) is stable against a variation of the Borel parameter T. The 
latter is varied in the range where both the hadronic and QCD representations of the 
correlator can be computed reliably. On the QCD side of the sum rules the restriction 
on T is mainly due to the perturbative contribution because T is an effective scale of 
as in Eq. (p^ . Taking into account the large value of the second order nonlogarithmic 
coefficient in Eq. (|2D|) we conclude that T cannot be chosen much less then 1.5 GeV where 
as(T)/7r ^0.1 to ensure the convergence of the perturbative series. Note that the power 
suppressed terms become dangerous at essentially lower T and that the above restriction 
also provides the convergence of the OPE. On the opposite side of the sum rules the 
hadronic representation of the correlator is reliable only for T ^ ujc which provides the 
exponential suppression of the contributions from higher resonances. 

The logarithmic dependence of npt(2;T) on T in Eq. ( plD is quite important for the 
stability of the sum rules. Therefore it is crucial to use the HQET renormalization group 
to get control over the high order logarithmic contributions. The leading and next-to- 
leading logarithms of T can be summed up by setting = T in the factor X{fi) and 
in the correlator npt(-2T). We adopt this prescription in our analysis. However, the 
normalization scale of in the 0{a'^) part of Ilpt{zT) is not fixed in our approximation 
and the corresponding /i-dependence is not compensated by X{fi). We do not use T as the 
normalization scale here when determining the optimal value of Uc because the resulting 
spurious T-dependence leads to rather unstable sum rules. If the normalization scale of a.^ 
in the O(a^) contribution is not correlated to T the result has a rather weak dependence 
on /i when varying the latter in the same interval as T. 

Adopting the central values of the input parameters and fi = 2 GeV, we obtain for 
the universal HQET decay constant 

fp = 410 (MeV)^/^ , (40) 

which constitutes an average for 2.05 GeV < cDc < 2.10 GeV. For these values the highest 
stability is observed. This value should be compared with the result obtained using the 
order expression of the correlator Hpt(tl') which is fp = 418 MeV at the optimal 
value Uc = 2.4 GeV. One notices that the inclusion of the 0{a'^) contribution leads to a 
rather small variation of fp though the correction to the correlator itself is quite large (cf. 
Eq. (^)). This can be explained by a considerable compensation of the large corrections 
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T (GeV) 



Figure 1: The B meson decay constant J'b as a function of the Borel parameter T of the 
Laplace sum rules for different values of the threshold parameter (D^. From top to bottom 
the curves correspond to ojc 



2.35,2.3,2.25 and 2.2 GeV. 



to Upti^j) and rrib (which enters the analysis through A) and the change of the optimal 
value LJc- 



6.2 B meson decay constant 

Taking into account the mass suppressed contribution to the Laplace sum rules and per- 
forming the analysis along the line described in the previous section, we obtain for the B 
meson decay constant 

fB = 206 MeV, (41) 

which constitutes an average for 2.25 GeV < tDc < 2.3 GeV. To illustrate the stability of 
the sum rules with respect to the variation of the parameter T we plot in Fig. |1] as a 
function of T for various values of Uc- It can be seen that the curves for ujc ~ 2.25 GeV 
and oJc ~ 2.3 GeV provide the most stable results. Note that for these values of Uc the 
function fsiT) has a weakly pronounced minimum around T ?s 2 GeV. At the same time, 
the use of the two-loop 0{as) approximation of npt(il;) along with the NLO value of m?, 
gives Jb = 205 MeV at the optimal value Uc = 2.75 GeV. As in the heavy quark limit we 
observe that due to the compensation of the corrections the O^a"^) result is practically 
the same as the 0{as) one. This fact is a strong argument in favour of our treatment of 
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the bottom quark mass. The Ola"^) mass suppressed corrections, which are included in 
Eq. (|4l|) , reduce the value of fs by approximately 5 MeV. Taking into account the fact 
that the nonperturbative part of the QCD contribution is saturated by the leading quark 
condensate, which gives about 10% of the total QCD contribution, we conclude that our 
result is stable with respect to all the types of corrections to the sum rules. 

Note that no rigorous results concerning the high order behaviour of the perturbative 
series both for the T-resonance and 5-meson sum rules are available and the absence of 
sizable higher order perturbative corrections to the decay constant within our approach 
cannot be proven strictly even though it works well up to the Ola"^) approximation. 
However, our approach provides the convergence also in higher orders if the divergence of 
the perturbative series for the T-resonance and 5-meson sum rules is related to the use of 
the pole mass and is thus dominated by the renormalon contribution. This is because we 
effectively remove the pole mass from the analysis of the decay constant in favour of the 
experimental moments of the T-resonance spectral density and operate with the relation 
between physical observables which is free from the corresponding renormalon ambiguity. 
Let us next discuss the uncertainty of the result in Eq. (^T]). The error in rrib of 



±60 MeV ^ results in an uncertainty of ±12 MeV in fs- The variation of the 



input value of as{Mz) within the experimental error bars 0.1185 ± 0.0020 |3J] leads to 
the uncertainty interval 201 MeV < Jb < 213 MeV. On the other hand, the result is 
not sensitive to the normalization point of as in the contribution and to the non- 

logarithmic three-loop coefficient in Eq. (^). Note, that the change of the parameters 
requires every time a new optimization of Uc- 

Another source of errors is the intrinsic uncertainty of the method due to the approx- 
imation of the hadronic contribution to the dispersion integral (^. A rough estimate of 
this uncertainty is obtained by the variation of the upper bound of the duality interval 
around its optimal value. The variation of Uc by ±100 MeV leads to ±7 MeV variation 
of fs which can be read off Fig. ^ A larger deviation from the optimal value leads to the 
essentially unstable sum rules. If we add the errors induced by the uncertainties in nib, 
as and cDc discussed so far in quadrature we obtain = 206 ± 16 MeV. 

A more advanced way to estimate the intrinsic uncertainty of the sum-rule approach 
is to change the weight function in the dispersion integral and redo the analysis using the 
sum rules which operate with Hilbert moments of the correlator ( ]T0| ) instead of its Borel 
transform (j^). The range of n relevant for reliable predictions of Hilbert sum rules is, in 
fact, quite restricted. The contribution of the mixed condensate grows rapidly with n. 
Thus, requiring the convergence of the OPE sets an upper limit on n. For the bottom 
quark it is n ^ 12 where the contribution of the mixed condensate is approximately 2/3 
of the leading quark one. At the same time, to avoid strong dependence of the result on 
ujc, one should use n > rrth/Cjc ^ 3. The optimal value of Cjc can be found by minimizing 
the dependence of the result on n in the above interval which results to Uc ~ 2.2 GeV. 
For this value the decay constant stays within the interval 192 MeV <, fs ^ 195 MeV as 
n varies from 4 to 12 as can be seen in Fig. |^ where fs is plotted as a function of n. This 
result for fs is in good agreement with the value obtained from the Laplace sum rules. 

Note, that in addition to the stability requirement there is a strong consistency check 
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Figure 2: The B meson decay constant fs as a function of the moment number n of 
the Hilbert sum rules for different values of the threshold parameter Qc- Ftom top to 
bottom the curves correspond to cDc = 2.35, 2.3, 2.25, 2.2, 2.15, 2.1 and 2.05 GeV. The 
value fjL — 2 GeV has been adopted. 

on the value of ujc- The decay constant drops out in the ratio of two moments which leads 
to a relation between Mb, rrib and Cuc- Our analysis is consistent if the physical value 
of the meson mass is reproduced from this relation for some value of the bottom quark 
mass in the interval given by the sum rules for the T resonances. For the above value of 
Uc this requirement is fulfilled for all n in the allowed interval. In fact, the ratio of the 
6*^^ and 7*^^ moment and the ratio of the 10**^ and ll*'^ moment imply exactly the central 
value nil, = 4.79. Note that for these moments fs has a local cxtremum as a function of 
n. In other words, if the 5-meson sum rules in the three-loop approximation are used 
to determine m^, the result is in perfect agreement with the NNLO value of the bottom 
quark mass from the T-resonances sum rules. 

We would like to emphasize that the dependence of Jb on rub is completely differ- 
ent for the Laplace and Hilbert sum rules. Thus, by comparing the results of the sum 
rules one can also estimate the error due to the uncertainty in m^. Furthermore, the 
prescription how the parameter ujc is determined and the structure of the condensate 
contributions are also quite different. The fact that both approaches give close results 
reflects the small intrinsic uncertainty of the sum-rule method in this particular case and 
furthermore strongly supports our treatment of mf,. Thus the total error originating from 
the uncertainty in rrib and the approximation of the hadronic spectrum can be estimated 
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as ±15 MeV. The remaining error is mainly due to the uncertainty in the input values of 
as{Mz). Consequently, as a conservative estimate of the uncertainty of our result for fs 
we quote ±20 MeV. Because of the stronger dependence of the Hilbert sum rules on rrib 
we use them to estimate the error but take the central value for our final result for Jb 
from the Laplace sum rules. 



6.3 D meson decay constant 

Since rric is not large in comparison to the scale A the expansion in l/rric can not provide 
us with the same accuracy as we have for the bottom quark. By the same reason no 
resummation of the charm quark mass logarithms is necessary and we can just use the 
three-loop result for the correlator up to in the full theory |T9|, 0. The use of the 

Laplace sum rules then leads to 

Jd = 195 MeV, (42) 

for the optimal value lJc = 2.35 GeV. The Hilbert sum rules are not reliable in this 
case because of the strong dependence on rric which is known with much less relative 
accuracy than nib- Furthermore, the Hilbert sum rules suffer from large contribution 
from the higher dimension condensates. The variation of the input parameters basically 
leads to similar variations of /d and Jb- Due to the weaker sensitivity of to the charm 
quark mass the additional uncertainty in does not lead to a larger error in Jb,. Thus, 
assuming the same intrinsic uncertainty of the sum rules we obtain the same error bars 
for the extracted value of fo- 



7 Conclusions 

To summarize, we have computed the B and D meson decay constants within the QCD 
sum rules approach. Our final results read 

/b = 206 ± 20 MeV , (43) 
fo = 195 ±20 MeV. (44) 

For the analysis we used the tree-loop result for the heavy-light current correlator. The 
large logarithms of the bottom quark mass have been taken into account by means of 
the HQET renormalization group. The bottom quark mass which essentially limits the 
accuracy of the sum rules for J'b has been extracted from the T-resonance sum rules up 
to NNLO. In the case of the B meson the analysis has been performed by employing 
Laplace and Hilbert sum rules. They have quite a different structure especially as regards 
the dependence on the bottom quark mass. The fact that the results obtained with these 
two approaches are in a good agreement gives us confidence in the reliability of the sum- 
rule method apphed to the calculation of the decay constants. The result also turned out 
to be quite stable with respect to inclusion of the perturbative corrections in as and l/m^ 
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and the nonperturbative corrections due to the vacuum condensate contributions. This 
allowed us to reduce the uncertainty of the extracted value of fs and fr,. The obtained 
values of both the B and D meson decay constants are consistent with the existing sum- 
rule results (cf. Eqs. (||) and (|D). However, the accuracy of our result is increased in 



comparison to the previous estimates. The values in Eqs. ( ^ ) and ( PD are in impressive 
agreement with the results obtained in lattice calculations (cf. Eqs. (|l|) and @). Probably 
no further improvement of the accuracy is possible within the standard QCD sum-rule 
framework due to the intrinsic uncertainty of the method. 

Our final comment concerns the current experimental status. A measurement is only 



available for the Dg meson decay constant where the most recent result reads [^] 

f^P = 286±44(stat) ±41(syst) MeV. (45) 
Converting the number in Eq. ( PD with the help of the lattice result fos/fo = 1-18 p8| , 



which agrees with the values given in |17|, we obtain fo^ = 230 MeV. This is in reasonable 



agreement with the experimental value of Eq. 
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